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ABSTRACT 

A ship berthing plan assigns surface vessels a berth 
prior to their port entrance, or reassigns ships once in port 
to allow them to accomplish in a timely manner maintenance, 
training, and certification events which build readiness for 
future operational commitments. Each ship requires different 
services when in port, such as shore power, crane services, 
ordnance, and fuel. Unfortunately, not all services are 
offered at all piers. At present, ship berthing plans are 
manually prepared by a port operations scheduler and often 
result in unnecessary berth shifts, which puts ships out of 
action for several hours. 

An extensive user-friendly computerized optimization 
model is developed and tested to assist the schedulers in the 
creation of a berthing plan which minimizes port loading 
conflicts, thus promoting fleet readiness through berthing 
stability. Norfolk Naval Station is used as an example 
because it exhibits all the richness of berthing problems the 
Navy faces. 
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I . INTRODUCTION 

While most ships' missions are executed at sea, their 
inport time is essential to maintain a high degree of 
material readiness and crew morale. A key to maintaining 
this is an efficient ship berthing plan. A ship berthing 
plan assigns surface vessels a berth prior to entering port 
or reassigns ships once in port "to accomplish a progression 
of maintenance training and certification events which build 
readiness for future operational commitments". (Wing, 1986, 
p.8) These events include, but are not limited to, combat 
systems maintenance, tests, and training, amphibious inport 
deck evolutions and other inport functions relevant to an 
individual ship class (COMNAVSURFLANT, 1987). In this study, 
the Navy's largest Naval Base is analyzed and modeled: Naval 
Station Norfolk, Virginia (NAVSTA NORVA) , exhibits all 
features required by other bases. A computerized 
optimization model is developed and tested to assist the 
schedulers in the creation of a berthing plan which minimizes 
port loading conflicts, thus promoting fleet readiness 
through berthing stability. 

The mission of the U.S. Navy, as set forth in Title 10, 
U.S. Code is: 

. . .to be prepared to conduct prompt and sustained combat 
operations at sea in support of U.S. national interests; 
in effect, to assure continued maritime superiority for 
the United States. (NWP-1, 1978. p. 1-3-1) 
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vice Admiral John S. McCain, Jr., U.S. Navy, emphasizes 
in •'Command at Sea", that the Commanding Officer of a United 
States naval ship must understand his ship's particular goal, 
mission, personnel and readiness as well as the Navy's over- 
all objective and mission. (Cope, 1966, p.vii) 

"Supporting military strategy involves ... having units 
properly manned, trained, equipped, and supported." It is 
the shore establishment's responsibility to "support the 
operating forces in terms of personnel, material, supply, and 
fiscal procurement; training; maintenance; and planning and 
operational guidance." (NWP-1, 1978, pg. 11-2-1) 

An inport period achieves several goals (Wing, 1986, 
P-8) : 

1. Enhance material condition of the ships through periods 
of maintenance in port at the unit (shipboard) , 
intermediate, and shipyard levels; 

2. Ensure crew proficiency through formal shore-based 
training ; 

3. Certification of public and crew safety and crew 
proficiency in the operation of installed equipment 
and systems ; 

4 . Provide adequate homeport time between operational 
periods in order to enhance morale; 

5. Conduct inspections and certifications mandated by 
public law. 

Prior to the port arrival of a commissioned naval ship or 
fleet auxiliary ship, the Commanding Officer sends a message 
to the appropriate naval authority stating the logistic 
requirements (LOGREQ) of his ship during the period in port 
(NWP-7, 1983, p. 7-1). This LOGREQ specifies any requests a 
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ship may have due to upcoming inspections, operational 
commitments, maintenance requirements or any other 
consideration the Commanding Officer identifies. 

Port operation ship berthing schedulers review logistic 
recjuirements , quarterly employment schedules and squadron 
requests for all home-based and visiting ships, and make 
berth assignments based on fleet requirements and port 
capabilities. Factors considered in berth assignments 
include: pier service requirements, deployment status, 
special operational tests, ship and berth characteristics, as 
well as crane requirements for on- or off-loading supplies. 

These considerations must be taken into account since 
each berth is unique in its capabilities: for instance, 
shore power and crane services available, depth and length of 
slip, fuel or ammunition loading capability and fendering 
system. (Papworth, 1988) 

An optimal ship berthing plan which minimizes port 
loading problems would require that all possible berths for 
each vessel be examined and "the one which best promotes 
fleet readiness while minimizing conflict between the inport 
goals would be chosen." (Wing, 1986, p.9) As a practical 
matter, this is impossible for a human scheduler to do. 
There are simply too many details to consider over time, and 
comparison of the "goodness" of alternate plans is 
problematic. 
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In order to assign ships to berths that offer required 
services while minimizing the number of berth shifts 
required, a high-speed computer should be utilized, berthing 
rules and ship priorities formally quantified, and an 
appropriate measure of effectiveness developed. The model 
developed herein meets these criteria. (Wing, 1986, p.9) 



A. PROBLEM SCOPE 

The focus of this thesis is on berthing surface ships 
assigned an inport period at the Naval Station Norfolk, 
although the methods developed here may be extended to other 
bases and stations. 

The mission of Naval Station Norfolk is 

...to provide, as appropriate, logistic support for the 
Operating Forces of the Navy, and for dependent 
activities and other commands as assigned. . . .The Port 
Services Officer (also known as the Port Operations 
Officer) is responsible to the Naval Station Commanding 
Officer for the performance of the port services 
functions. ...For ships (units) under naval control, the 
port service's functions include the assignment of berths 
and anchorages; the use of piers, landing sites, pilots; 
coordination of logistic requests for supplies, fuel, 
medical services, communications, hazardous material 
handling and other services.... (Fleet Guide, 1986, p. 5- 
3) 



The Norfolk Naval Station consists of 15 piers, depicted 
in Figure 1, which exhibit different physical characteristics 
and services. Typically, the average number of ships in port 
is 50 with the highest port load peaking at 74 during the 
Christmas holiday. These vessels usually rely on shore power 
rather than on their own power. (Papworth, 1988) Shore power 
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and other facilities permit ships to operate and test combat 
systems and other mission capabilities while in port 
( COMNAVSURFLANT , 1987) . The increasing number of ships 
homeported at Norfolk (presently 118) , along with unique 
requirements by ships and lengthy pier maintenance projects, 
combine to make pier scheduling an extremely difficult task 
requiring complex planning (COMNAVSTANORVA, 1987) . 
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Figure 1. Naval Station Norfolk Piers 

Piers have different lengths, water 
depths, crane services, and power, etc. 
available. 
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B. CURRENT PROCEDURES 



The Naval Station Norfolk ship berthing plans are 
manually prepared by the schedulers with the aid of pen and 
paper and a wall-size mock-up of the pier layout with scale- 
size ship silhouettes. Once informed of which ships are 

scheduled to be in port for the next week, the scheduler 
first determines which berths can physically accommodate each 
ship. 

The berth scheduling rationale is based on the following 
primary criteria: 

1. The ship's length must be less than the length of the 
pier 7 

2. The pier-side depth must be five feet greater than the 
ship's draft to allow for tidal change as well as 
propeller wash and engineering plant requirements; 

3. The ship's beam plus fender system must extend less 
than the distance between the berth and the next 
closest pier or berthed ship plus room to allow a ship 
to maneuver; 

4 . The berth should provide at least the minimum required 
number of shore power cables (COMNAVSURFLANT, 1987) . 

After the physically feasible berths are identified for 
each ship, the scheduler then considers a secondary set of 
guidelines specific to Norfolk listed in Appendix A. Every 
port has an analogous set of local berthing criteria. 

At this point, scheduling becomes difficult. The 
scheduler assigns each ship to a feasible berth and tries to 
maximize the proportion of requested services and minimize 
the number of berth shifts required to accommodate these 
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requests over time. This berth plan is the initial input to 
a key planning event, the berthing conference. 

Once a week, a berthing conference is held at port 
operations and attended by representatives from squadrons, 
groups, type commanders. Military Sealift Command, Norfolk 
Supply Center, Public Works Center (PWC, utilities and crane 
scheduler) , Readiness Support Group and Port Operations 
(scheduler, chief pilot, ordnance officer, dockmaster and 
policy maker). The squadrons all represent their ships' 
requests for docking and undocking times, as well as for 
particular berth assignments. PWC advertises feasible pier 
utility services. The pilot assigns move times for ships 
constrained by tide. Compromises are worked out and the Port 
Operations Officer makes final decisions. (Papworth, 1988) 

The final berthing plan resulting from the berthing 
conference is used as the start of the following week's 
schedule. Coordination among all these participants is 
vital. Changes in the announced plan are inevitable — the 
schedule often changes hourly. The sheer frequency of 
revisions makes a strong case for the use of a computerized, 
optimizing berthing plan. The consequence of oversights is 
delay, and delays cost time and money. 



8 



II. AN OPTIMIZING BERTHING MODEL 



The goal is to create an optimal berth plan, at a daily 
level of detail, for all ships scheduled to be in port during 
the prospective planning horizon (say, a week) . As 
demonstrated in Chapter I, the berthing requirements are well 
defined. This chapter explains the basic model developed to 
satisfy these requirements and produce optimal berth 
schedules . 

The situation calls for a set of discrete ship-to-berth 
assignments, with limitations on feasible assignments. These 
limitations (on length, draft, power cables, and so forth) 
are easily expressed as linear functions of ship-to-berth 
assignment variables. This suggests a linear integer 
program. 



A. LINEAR INTEGER PROGRAM MODEL FORMULATION 



Indices : 
i = 1 , . . . , I 

j = !,...,! 
P = 1, . . . ,P 
b = 1 , . . . , B 
n = 1, . . . ,N 

t = 0 , . . . , T 
k = 1, . . . ,K 
q — 1,...,Q 



individual ships 

individual ships (alternate index) 

piers 

berths 

nesting position (l=pierside, . . . ,N-1= 
next-to-last position) 

day (t=0 indicates current day) 

basin 

services 



9 



m = 1 , 


. ..,M ship characteristics (draft, 

length, . . . ) 


z = 1 , 


. ..,Z pier characteristics (depth, 

length, . . . ) 



For simplicity of presentation, it is implicitly 
understood in the following that only permissible 
combinations of indices are used. 

Given and Derived Data: 



^pk 


1 if pier p belongs to basin k; 0 otherwise 


L 

P 


length of pier p + allowable overextension of 
ships during high port loading 


1 . 
1 


length of ship i + minimum distance between 
adjacent ships (bow/stern, stern/stern, 
seawall/ ship) 


E 

P 


number of power cables available at pier p 


®i 


minimum number of power cables required by 
ship i 




width of basin k - tug maneuvering room 


'^i 


beam of ship i + fendering 


®i 


N if ship i cannot physically nest; 
1 otherwise 


^i 


1 if ship i can berth ships outboard; 
0 otherwise 


"ip 


1 if fendering and superstructure on pier p is 
compatible with ship i; 0 otherwise 


SD. 

im 


characteristic m for ship i 


PD 

pz 


characteristic z for pier p 


R 

n 


reward for nest position 


LQ. 

^iq 


ship i priority for requested service q 
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1 if Pier service q is available on pier p; 

0 otherTtfise 

safety distance between ship draft 
and water depth 

1 if t> 1 and ship i can physically fit in 
specified berth b at pier p, at nesting 
position n, and is scheduled to be in port on 
day t; 

0 otherwise, in particular when s. = N and 
n > 2 

to help the human scheduler, rather than 

(foolishly) try to replace him, extensive user-friendly 

facilities should be provided to allow manual assignment of a 

ship to a specific berth, subset of piers/berths, or nesting 

position. These coercions are simulated in the prototypic 

implementation via input of derived compatibility data, 

C. ... This allows the scheduler to restrict any or all 
ipbnt 

permitted indices for a ship, i.e., a specific berth, group 
of berths/piers, and/or nesting position for a specific ship 
during any or all days the vessel is scheduled to be in port. 
When the user identifies specific requests, all other 
are automatically set to zero, thus ensuring the ship will be 
berthed only as specified by the scheduler. 

A ship may be assigned to one of the specified berths at 
a pier as long as all of the primary berthing conditions (1)- 
(5) are satisfied. If these primary berthing criteria are 
violated for every pier associated with each specified berth, 
the ship can not berth and the problem is deemed infeasible. 



qp 

DR 



^ipbnt 



In order 



11 



Ur aft ^ P depth ~ DR 



(I) 




( 2 ) 



^ D i arrive ^ ^ 



I, arrive 



(3) 




( 4 ) 




( 5 ) 



Condition (1) ensures the pier depth is deep enough for 
the ship's draft plus safety distance. Condition (2) berths 
a ship only if its length does not extend past the pier. For 
a ship to be considered compatible, it must be scheduled to 
be in port during the day considered as ensured by conditions 
(3) and (4). Condition (5) does not allow a ship to be 
assigned a berth where it would have a fendering or 
superstructure interference. 

The objective is the "goodness" of any given feasible 
berthing plan. The problem is greatly simplified if this 
benefit can be expressed as an additive, separable linear 
function of individual ship-to-berth assignments. To provide 
such an objective function, individual ship service requests 
are prioritized among and between ship classes: larger ships 
such as aircraft carriers are authorized higher priorities 
for services than destroyers or frigates. The benefit is 
expressed as a function based on this ship priority for 
services, pending inspections, deployments, whether the pier 
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offers any or all of the requested services and how far into 
the future the decision will be committed. 



Recognizing the time value of information and 
uncertainty, an exponential function discounts the preference 
awarded to a ship desiring a berth in the future vice a ship 
requesting it today. 



The benefit of a potential assignment is thus calculated 
by summing term-by-term the pairwise products of the weighted 
ship requests (LQ) with the vector (A) which identifies 
services available at each pier. This is an indication of 
how well each berth satisfies a ship's needs. The inspection 
and deployment (SD) factors are then added to the weighted 
ship requests; this allows a ship with an upcoming inspection 
or deployment to be ranked higher than other ships of the 
same type. The updated weight is multiplied by an 
exponential term to give a greater consideration to ships 
requesting services today than ships scheduled to be in port 
in the future. Lastly, a reward (R) based on nesting 
position is added to yield the final benefit for assigning 
the ship to a specific berth. This final nesting position 
reward encourages the model to berth ships pierside. 



BN. 



ipbnt 



benefit from berthing ship i, at pier p, in 
berth b, at nesting position n, on day t; 
derived only if C, , . = 1, and defined as 

follows: ^ 



ipbnt 
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Variables : 



X. . . A binary variable specifying if ship i is to 

^ be berthed at pier p, in berth b, at nesting 

position n, on day t. In the implementation, 
the variable, X. ^ , is included in the model 
only when the corresponding assignment 
satisfies the feasibility conditions 1-5. 

To take into account the fact that the 
berthing of ships is an ongoing process, t=0 
indicates a ship's current position. Figure 2 
illustrates the meaning of each index 
graphically. 




Zi bnt indicates if ship i shifted to pier p, in 

^ berth b, at nesting position n, on day t. 

This variable is generated only if the ship 
was berthed on day t-1. 
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Formulation: 



Maximize ^ ipbm^ipbnt X penally^ ipbnt 

ipbnt ipbni 

Subject to: 



Ji^ipblt ^ ^p p,l e (PxT) (6) 

lb 

ipbnt ^ ^p P,t e {PxT) (7) 

ibn 

^h^ipbnt ^ X^'"V(«+i)r e {PxBx{N -\)xT) (8) 

/ / 

pbn 

'Y^^ipbm ^ i e {PxBxNxT) (10) 

i 

ipbnt {PxBxT) (11) 

in 

{N - fl)Xipbnt + '^Xipb[n+l)t + '^^ipb(n+2)t ^ X - fl ( 12 ) 

j^i J¥=i 



i,p,b,n,t 6 (/|(u, = \)xPxBx{N — \ )xT) 
'Z»iDptX,pt„ ^ k,b,te(KxBxTe(SDp„,p.,<l)) (13) 

ipn 

^ ipbnt ~ ^ipbn(t—\) ^ipbnt ^ ^ 
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i,p,b,n,t e {IxPxBxNxT) 



(14) 



^ipbnt ^ {^) ^ } 



i,p,b,n,t e. {IxPxBxNxT) 



( 15 ) 



^ipbrtt'^^ i,P,by(i,t e (IxPxBxNxT) ( 16 ) 

In the above formulation, the objective function is to 
maximize the ship-to-berth assignment benefits less a berth 
shift penalty. This penalty decreases the total benefit of 
the plan each time a vessel is required to move to a 
different berth or nesting position from day to day in order 
to receive required services at a new berth or to free its 
current berth for another ship. Since the formulation 
encompasses the entire planning horizon, the optimal plan 
takes into account the arrival on any day of new ships and 
their required services. Initial ship positions are treated 
as arrivals on day 0. 

Constraints (6) require that the total length of ships 
berthed at pier p and nested inboard are less than the length 
of the pier plus allowable extension. Constraints (7) ensure 
that each pier has sufficient number of power cables to 
support ships berthed alongside. Constraints (8) allow ships 
to be nested outboard another vessel only if its length is 
less than or equal to that of the inboard vessel. 
Constraints (9) ensure each ship is berthed at only one slip 
when scheduled to be in port while constraints (10) allow 
only one ship per berthing position. Constraints (11) ensure 
specific ships do not nest and constraints (12) preclude 
berthing outboard of incompatible ships. Constraints (13) 
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provide room for a tug to maneuver among berthed ships in 
each basin. Berth shifting is calculated with constraints 
(14). Conditions (15) ensure the assignment variable is 

binary while (16) requires the berth shift variable to be 
nonnegative. 
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III. SHIP BERTHING PLAN EXAMPLE PROBLEM 



A prototype model has been evaluated using a GAMS 
generator (Brooke, Kendrick, and Meeraus, 1988) and MPSX 
solver (IBM, 1972) . The model has been tested using an 
example with seventeen ships, eight piers (see Figure 3), and 
a six-day planning horizon. 




Figure 3 . Norfolk Naval Station 

Piers 4S, 7, 10, 11, and 12 are scheduled in the 

example problem. 

This example incorporates a wide variety of ship types: 
carrier, frigate, destroyer, cruiser, oiler and battleship. 
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The physical characteristics input for each ship include: 
length, draft, beam, arrival date, departure date, number of 
shore power (PWR) cables required, whether the ship can nest 
(SHP) or allow ships outboard (NOOUT) . Upcoming inspections 
or deployments are identified along with the penalty incurred 
if a berth shift is required (COSTSHF) . Table 1 displays a 
sample of the GAMS data input for the seventeen ship example 
problem. 



TABLE 1. 

EXAMPLE PROBLEM SHIP CHARACTERISTICS 



SHIPDATA (I , 


pSDATA) 






LIST OF 


SHIPS AND CHARACTERISTICS 






LENGTH DRAFT 


BEAM , 


ARRIVE 


DEPART 


SHIPPWR 


INSPECT 


DEPLOY 


SUP 


COSTSHF 


AFS2 


581 


24 


79 


5 


5 


4 








400 


A0R4 


659 


33.3 


96 


3 


5 


3 








500 


ODG-I 


437 


20 


47 


2 


5 


3 








300 


DOG 6 


437 


20 


47 


1 


5 


3 








300 


LIIA-t 


840 


26 


106 


2 


5 


14 






1 


400 


enei 


887 


38 


1 08 


1 


2 


6 








500 


Bnr>ia 


887 


38 


108 


5 


5 


6 








500 


CVN71 


1300 


37 


134 


1 


5 


8 






1 


1000 


CV67 


1300 


35.9 


130 


1 


5 


24 




200 


1 


1000 


CV66 


1300 


37 


130 


1 


3 


24 




200 


1 


1000 


CG2 7 


547 


28.8 


54.8 


1 


4 


4 








350 


CG3 0 


547 


28 . 8 


54 .8 


2 


5 


4 








350 


CG34 


547 


28.8 


54 . 8 


4 


5 


4 








350 


CG4 8 


567 


33 


55 


1 


3 


6 


100 






350 


CG51 


566 


31 


55 


1 


2 


6 








350 


FFG5 


414 


24.2 


44.2 


1 


4 


2 




100 




300 


TAF8 


524 


22 


72 


1 


2 


2 








400 


TA0189 


678 


34.5 


97.5 


2 


4 


4 








4 00? 



To identify any particular ship or ship type, refer to 
Jane's Fighting Ships . 1988. Each pier is characterized in 
Table 2 by its length, depth and shore power available: each 
basin width is identified along with the piers that form the 
basin . 
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TABLE 2. 

PIER CHARACTERISTICS 



PIERDATA(P, PDATA) LIST OF PIERS AND CHARACTERISTICS 





riERL 


DEPTH 


POWER 


3 2N 


1300 


50 


24 


12S 


1300 


50 


24 


13N 


1397 


50 


24 


US 


1397 


50 


24 


ION 


1300 


38 


56 


7N 


1350 


45 


24 


7S 


1350 


45 


21 


4S 


1347 


40 


2 4 ; 



The services available pierside include: diesel fuel 
(DEM), JP5 fuel, Military Sealift Command (MSC) , Stores, 140T 
crane, DESR0N2 (DRON2) and C0MDESGRU8 (CDG8) sponsorship, and 
ordnance handling certification. Table 3 shows the pier and 
service availability GAMS input used in the sample problem. 



TABLE 3. 

EXAMPLE PROBLEM PIER/SERVICE AVAILABILITY 



AVAIL(P,J) 



LIST OF PIERS WITH SERVICES OFFERED 



DFM JP5 MSC STR 140T 4160V DROH2 ORDN CD68 TEND 



32N 1 1 1 

12S 11 1 

UN 1 

US 1 

ION 

7N 1 1 

7S 1 

4S 1 1 1 1 1; 



1 

1 

1 

1 

1 

1 



1 



The 

request 



maximum "weights" authorized 
services are listed in Table 4. 



for ship types to 
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TABLE 4. 

MAXIMUM WEIGHTS AUTHORIZED FOR SHIP/SERVICE REQUESTS 



SHIPS 


WEIGHTS 


CV/CVN 


1000 


BB 


900 


LHA/AOE/AOR/AD/AO 


800 


TAO/TAF/AFS 


700 


LPH/LPD/LST 


600 


CG/CGN 


500 


DD/DDG 


400 


FF/FFG 


300 


ARS/MSO 


200 



The maximum weight limits are used to edit the priority 
for services requested by a ship. 

The weighted values assigned to each ship for requested 
services are easily identified in the GAMS data matrix 
utilized in the example problem as seen in Table 5. 

TABLE 5. 

EXAMPLE LOGREQ PRIORITIES FOR SHIP SERVICES 



LOGP.EQ(I , J) 




WEIGHTED SHIP TO 


SERVICES REQUIREMENTS 




DFM JP5 


MSC SIR 


140T 4160V DRON2 ORDN 


CD68 TEND 


AFS2 




600 


600 




A0R4 


750 750 








LHA4 










DDG4 










DDG6 










BB61 










BB61a 


600 








CVN71 


900 




999 




CV67 


900 








CV66 


900 








CG27 








400 


CG3 0 








400 


CG3 4 








400 


CG4 8 








400 


CG51 








400 


FFG5 


200 








TAF8 




700 700 


400 




TA0189 


700; 
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For instance, BB61a indicates a second inport period for 
BB61 during the planning horizon; A0R4 has an assigned 
weight of 750 for the requested services of DFM and JP5 , 

The remaining physical characteristics for all ships, 
piers and basins essential to the problem are given in 
Appendix B. The resulting integer program is generated by 
GAMS and solved using the MPSX solver. The final GAMS output 
berthing plan is displayed in Table 6 and illustrated in 
Figures 4 through 9. Each daily berth plan is printed to 
show all ships scheduled to be in port and their assigned 
berth. The dashed silhouettes in Figures 4 through 9 

indicate a ship departure and the arrows identify berth 
shifts . 

TABLE 6. 

GAMS OUTPUT FOR FINAL BERTHING PLAN 



580 PARAMETER OAV • l.OOO COUMTtR OF DAY 



580 parameter sol 

12N.2.1N80 12S.1.1NBD nN.2.lM8D lOH.l.lMBD lOH.l.OUTBD 10N.2.1NBD 7H.1.1NBD 7S.I.INBD 7S.2.1NBD 4S.2.1NbO 



CVN7 1 

CV6 7 

CV66 

cg:7 

CG68 

CG51 

FFG5 

TAF8 

B86 1 

DDGt> 
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300 PARAMETER DAY 



r.OOO COUNTER OF DAY 



380 PARAMETER SDL 

IZN.r.lNBD 125.1.1N8D 11N.2.INBD llS.l.INBD 115.2. 1N8D 1 ON . 1 . INBD 1 ON ^ .DUTBD 1 ON .2 . 1N0D 1 ON. 2 . DUT8D ?N.1.INBD 



CVN71 
CVS 7 
CV6S 
cg:7 
CGJO 
CGC8 
CG51 
FFG5 
LHA« 
DDG< 



1.0 



1.0 



1.0 



1.0 



1.0 



1.0 



1.0 



1.0 



• 7S.I.1NBD 7S.2.1NBD AS.l.lNBD <S.2.1NBD 

TAFB 

TA01B9 1*0 

8BS1 1.0 

DDGS 1.0 



3B0 parameter day • 5.000 COUNTER DF DAY 



3B0 PARAMETER SOL 



12N.2.1NBD 12S.I.1NBD 11N.2.1N8D llS.l.INBD 11S.2.1N8D lON.l.lNBD lON.l.DUTBD 10N.2.1NBD 7H.I.1NBD 7N2.1NBD 



ADR< 

CVN71 

CVS 7 

CVSS 

CG2 7 

CG30 

CG<48 

PFCS 

LHAS 

DDG4 



1 .0 



I.O 



1.0 



1 .0 



♦ 7S.2.1NBD SS.l.lNBD 



TAD189 1.0 

DDGS 1.0 
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— 


5B0 TARAMETER DAV 


• 


4.000 COUNTER DP DAV 








— 


SBO PARAMETER SOL 














12H.2.1NBD UN.2.1N8D 


ns. 1. INBD 


ns. 2. INBD ION. 1 . INBD ION. I .OUTBD ION. 2. INBD 


7N. 1. INBD 


7N.2. INBD 


7S.2.1NBD 


A0R4 










1.0 




CVN71 


I.O 












cvt? 


1.0 












cg:7 






1.0 








CGSO 






1.0 








rrc5 








1.0 






LMA4 






1.0 








DDC4 




1.0 










DDG« 












1.0 


CG5« 






1.0 








• 


43. 1 . INBD 












TA0IB9 


1.0 













— 


SBO parameter day • 


S.OOO COUNTER DF DAY 












— 


SBO parameter SDL 
















12N.2.INBD 12S.2.1NBD nN.2.1HBD 


ns. 1. INBD ns. 2. INBD ION.] 


;. INBD 


ION. 2. INBD 


7N.2.1NBD 


7S.2. INBD 


4S.2.INBD 


APS2 














1.0 


AOR4 










1.0 






CVN71 


1.0 














CV67 


1.0 














CGiO 






1.0 










EB61A 


1.0 














LHA4 




1.0 












DDG4 




1.0 












DDG6 












1.0 




CGS4 








1.0 









Figure 4 shows the position of ships in port at start of 
the planning horizon. 
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Figure 5 illustrates that on Day 1 both FFG5 and CG27 are 
required to berth shift in order to make room for the arrival 
and berthing of CV67. The slashes outboard CG27 indicate the 
ship's request for none to berth outboard. 




CG51 berth shifts to allow CG30 to berth pierside in 
accordance with the scheduler's input. These activities are 
shown in Figure 6. 
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Figure 6. Day 2 Ship Berthing Sample Plan 
BB61, CG51, and TAF8 depart and are underway on Day 3 

whereas AOR4 arrives in port. When CG51 leaves its pier, 
CG3 0 berth shifts to be pierside vice nested out. These 
events are displayed in Figure 7. 




Figure 7. Day 3 Ship Berthing Sample Plan 



CG48 and CV66 are underway from Norfolk on Day 4. CG34 
arrives at the base and assumes CG48's berth illustrated in 
Figure 8. 
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Figure 8 . Day 4 Ship Berthing Sample Plan 



On Day 5 AFS2 arrives inport and BB61 returns but to a 
different berth. TA0189, FFG5 , and CG27 are underway for 
sea. Figure 9 identifies these movements. 




DAY 5 



Figure 9. Day 5 Ship Berthing Sample Plan 
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IV. COMPUTATIONAL EXPERIENCE 



GAMS, the General Algebraic Modeling System (Bisschop and 
Meeraus, 1982) , ”is designed to make the construction and 
solution of large and complex mathematical programming models 
more straightforward for programmers and more comprehensible 
to users of models.” (Brooke, Kendrick, Meeraus, 1988, p. 
xiii) GAMS has been developed to: 

1. Provide a high-level language for the compact 
representation of large and complex models; 

2. Allow changes to be made in model specifications simply 
and safely; 

3. Allow unambiguous statements of algebraic 
relationships ; 

4. Permit model descriptions that are independent of 
solution algorithms. (Brooke, Kendrick, Meeraus, 1988, 
p.3) 

Using GAMS to implement the prototypic ship berthing 
model enabled experimentation and easy changes to both the 
model and its data. However, when analyzing GAMS model 
results, the user must be very careful to keep the underlying 
mathematical model in mind: it is very easy to forget about 

mathematical programming theory and just concentrate on the 
rather hypnotic allure of powerful GAMS statements. This can 
(and did) confound verification (debugging) and validation. 

The example ship berthing plan problem was originally run 
using GAMS/MINOS (Murtagh and Saunders, 1983) to debug the 
model. Unfortunately, MINOS has no integer capabilities. 
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since the linear programming relaxation of the problem did 
not naturally solve with integer values for the binary 
variables, the FORTRAN-based mixed-integer programming solver 
XMP/ZOOM was introduced. (Marsten, 1981) and (Singhal, 
Marsten and Morin, 1987) 

The ZOOM solver is "intended for medium-sized problems 
with no special structure and up to about 200 zero/one 
variables." (Brooke, Kendrick, Meeraus, 1988, p. 225) Zoom 
occasionally obtains respectable integer solutions, but only 
after a long series of experiments adjusting the ZOOM 
parameters and options. ZOOM consumes enormous computer 
resources in these trial-and-error experiments. Many 
shortfalls of the ZOOM solver have been identified and 
referred back to the developer, Marsten, via Meeraus. ZOOM 
can not even solve some trivial test problems due to several 
apparent severe bugs and has therefore been deemed 
inefficient and unreliable. 

With the assistance of Professor Terry Harrison of the 
Pennsylvania State University, an electronic mail connection 
via IBM-BITNET has been established permitting the GAMS model 
to be transported to and solved on the IBM 3090-400 at PENN 
State using GAMS/MPSX (IBM, 1978) . GAMS/MPSX worked, 
requiring 70 IBM 3090-400 processing seconds and 7,143 
iterations to solve the example model with 1,779 constraints, 
2,864 variables, 1,512 binary variables and 16,453 non-zero 
coefficients . 
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Motivated by these experiences, Professor Harrison is 
preparing a GAMS/X-System (Brown and Graves, 1975) interface, 
which he and Brown will test on problems including such as 
that reported here. The goal is to show that the port 
scheduling model can be solved quickly and inexpensively at 
realistic scale on a modest computer (Bausch and Brown, 
1988) . This is important, because the full-scale Norfolk 
berth scheduling problem will require some advanced 
optimization techniques not yet available via GAMS. To see 
this, consider that with 24 piers, 144 berths, and 74 ships 
inport an average of five days over a seven-day berthing 
plan, up to 120,107 constraints and 53,280 binary variables 
may be required. 

The formidable size of the hypothetical port scheduling 
model problems is mitigated by numerous restrictions on 
permissible realistic combinations of indices (berths, ships, 
services, etc.). The dollar operator feature in GAMS 
"provides powerful and concise exception-handling 
capability." Explicit if-then-else statements constructed 
within an equation or assignment makes a program more 
manageable by decreasing the number of equations and 
variables generated. (Brooke, Kendrick, Meeraus, 1988, p.72) 

The entire GAMS/MPSX input listing for the example 
berthing problem is in Appendix B. 

GAMS is a powerful tool, but expensive to use in terms of 
computer resources. The example model requires 30 IBM 3090- 
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400 seconds just to generate the input for an optimizer. 
After solution, the simple report writing takes 4.5 seconds. 
By contrast, models of equivalent size and complexity are 
generated in a second, or so, by use of customized problem 
generators written in general-purpose compiled languages 
(e.g., FORTRAN). Such old-fashioned generators take longer 
to write and debug than GAMS, and are less easy to modify, 
but they generate with enormously improved efficiency. Given 
that the port schedulers will not likely have an IBM 3090 
super computer available soon, or be willing to wait hours 
for each solution, it seems likely that a more conventional, 
old-fashioned approach will be called for. 



31 



V. CONCLUSIONS 



Optimization-based berth scheduling is feasible and 
effective. The prototype introduced and developed here gives 
compelling evidence that a computer-based model can express 
the berthing problem concisely in easy-to-understand 
displays, and automatically produces berthing plans capturing 
an enormous amount of the realism and detail that make such 
scheduling a challenging manual chore. Better yet, the 
method developed here encourages human interaction. 

In the context of the proposed model, extensive user- 
friendly facilities can be accommodated to allow a port 
operations scheduler to manually assign a ship to a specific 
berth, subset of piers/berths or nesting position. The 
optimization model then completes the tedious details of the 
berth plan. Thus, the port operations scheduler can 
naturally express any "human judgement" issues and the 
optimization assures that high-quality berth plans are easily 
and quickly produced. 

This optimization program would also give the scheduler 
the flexibility to evaluate alternate "what if" berthing 
plans. In this role, quick-response identification of 
upcoming infeasibilities may be as useful as comparative 
evaluations of the relative merit of alternate plans. There 
is no current manual analog for this capability, nor is it 
likely that the manual time and effort will be available to 
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devote to much more than cursory analysis of schedule 
changes . 

This optimal berthing plan model could easily be adapted 
for other naval bases and include submarines. 

Bases and Stations Information System (BASIS) is a new 
data base management system currently being installed for 
prototype tests at Naval Station, Norfolk. BASIS is planned 
for future use by all U.S. Naval Stations and Bases world- 
wide. BASIS is "a computer-based network developed to 
fulfill the command staff's need for timely and accurate 
information concerning base/station activities." (NPRDC, 
1988, p. vii) 

The Port Services module of BASIS is designed to support 
management in the complete and efficient supervision of Port 
Services functions by maintaining data that can be used in 
assessing current and future Port Services needs. The system 
"operates in an on-line environment utilizing Video Display 
Terminals that allow for interactive processing of data via 
add, change, delete and display functions." (NARDAC, 1988, 
p. iii) 

The Port Services module is capable of handling 
waterfront functions such as current and projected berthing 
plans, maintenance schedules, waterfront status, ship and 
pier characteristic data, ship schedules, pier maintenance 
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projects, and a variety of statistical data. (NARDAC, 1988, 
p. iii) 

An optimization-based berth scheduling module can be 
embedded in BASIS, providing a powerful decision aid to port 
operations management. Given the design specifications of 
BASIS, an optimal berthing model should be able to retrieve 
the information it needs. Following optimization, the video 
display features of BASIS would be invaluable in expressing 
the current and planned port activities. The pictorial 
display of a ship berthing plan is very useful and easier for 
the scheduler to interpret than an assignment table. Thus, 
the port operations scheduler could produce efficient and 
current berthing plans. 

Port scheduling is crucial to the U.S. Navy. Considering 
the tempo of schedule changes and the meticulous detail which 
preparation of every schedule must consider, a manual 
scheduler is hard-pressed to weigh myriad alternatives and 
fine-tune every alteration. It is inevitable that oversights 
will lead to delays. If an automatic, optimization-based 
decision support system prevents unnecessary delays or berth 
shifts, then such a system clearly contributes to the 
readiness of the fleet. 
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APPENDIX A 



NAVAL STATION NORFOLK BERTH SCHEDULING GUIDELINES 



1. Due to pier superstructure, the following ship types 
can not berth at these prohibited locations: 



SHIP TYPE 

LSD,LPD 

CV/CVN 



PROHIBITED LOCATIONS 
PIER 4 berth 5 and 6 
PIER 2,3,4,10 



2. The fendering system limits the ship types certain 
piers are or can be configured for. All other ship types may 
go to any berth provided it is physically feasible and shore 
power is available. 



SHIP TYPE 
BB 
LHA 
LPH 
LPD 

CV/CVN 



COMPATIBLE PIERS 
all except 10 
5,7,10,11,12,25N 
2,5,7,11,12 
2,3,4,5,7,10-5,11,12 
7N, 11,12 



3. Ships would like to be berthed at piers that their 
respective squadrons "sponsor". 



PIER 

20 

21 

24.25 

10.25 



SPONSOR SQUADRON 
SERVGRU4 
DESRON 2 
DESRON 10 
CRUDESGRU 8 



4. During high port loading, ships berthed bow out may 
extend up to 20 feet beyond the end of the pier. 

5. Maintain a distance of 50 feet between ships berthed 
bow-to-bow, bow-to-stern, stern-to-bow, and 25-50 feet 
between a ship's bow-to-stern and a seawall. 

6. The larger the ship, the higher its priority should 
be in receiving services. 

7. Do NOT nest CV, CVN, LHA, LPH and LPDs due to their 
hull structure. 



8. An outboard ship's length must be less than or equal 
to the inboard ship's length. This minimizes the stress on 
mooring lines. However, during high port loading, the 
outboard ship may be up to 20 feet longer. 

9. Preferably, berth ships in "UPKEEP" near a tender or 
SHip InteiTTiediate Maintenance Facility (SIMA) , responsible 
for repairs. 
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10. Certain services are rendered only at specific piers: 
e.g., refueling pier side, ordnance transfer, major stores 
loading, collimation (piers 5, 7, 24, 25; berths 1 and 2), 
sonar testing (bow out, end of pier), and cranes. 

11. The maximum number of ships nested is usually two but 
may go up to four. This is primarily due to shore power 
1 imitations . 

12. If LHAs require lowering their ramp, they must do so 
on piers 5,7,11 or 12. (The drive-on and -off ramp is used 
to load vehicles.) 

13. Certain ships MUST go to specific berths, (e.g., USS 
Mount V'Jhitney, Pier 25-1) 

14. Ships preparing for deployment and inspections have a 
higher priority for services than others of the same ship 
type. 

15. Two ships of certain classes may berth Chinese (bow- 
to-stern) . (Spruance, Oiler, BB, Ticonderoga, Yorktown, DDG, 
FF) This is not a major factor but may be a consideration. 
(This is an infrequent event.) 

16. Ships undergoing a Radiation Hazard (RADHAZ) survey 
must be 200 feet out of range of any line-of-sight shore 
structure or other ship's superstructure. (This is an 
infrequent event.) 

17. Berth ships (AOE, AO, AOR, AFS) requesting inport 
underway stream qualification training and tests, underway 
replenishment standard qualification trials (UNREP SQTS) 
across from each other in the same basin or across an 
unobstructed pier. (This is an infrequent event.) See Figure 
4 . 




Figure. 4 UNREP SQTTS Berthing Positions 
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18. SUBRON 6 and 8 are responsible for assigning 
submarines and tenders to berths on piers 22 and 23. Thus, 
piers 22 and 23 are removed from our considerations. 

19. Phone lines, fresh water, 125psi steam, and 
collection, holding and transfer (CHT) hook-ups are not scare 
at the Norfolk Base and are therefore not considered in the 
model . 
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APPENDIX B 



GAMS INPUT FOR SAMPLE BERTHING PROBLEM 

$TITLE SHIP SHIP BERTHING - LINEAR PROGRAMMING 

$OFFUPPER 

$OFFSYMXREF 

$OFFSYMLIST 

$OFFUELXREF 

$ONTEXT 

The model now will assign a ship to a berth during the time period 
the ship is scheduled to be inport. 

$OFFTEXT 



* INDICES 



SETS 

I SHIPS 

/AFS 2 , AOR4 , CVN7 1 , CV67 , CV66 , CG27 , CG3 0 , CG4 8 , CG5 1 , FFG5 , TAF8 , 
TA0189, BB61, BB61a , LHA4 , DDG4 , DDG6 , CG34/ 

SDATA SHIP DATA 

/LENGTH, DRAFT, BEAM, ARRIVE, DEPART , SHI PPWR, INSPECT, 

DEPLOY , SHP , COSTSHF/ 

P PIERS 

/12N, 12S, IIN, IIS, ION, 7N, 7S, 4S/ 

PDATA PIER DATA 

/PI ERL, DEPTH, POWER/ 

J SERVICE 

/DFM, JP5,MSC,STR, 140T , 4 160V , DRON2 , ORDN , CD68 , TEND/ 

T TIME PERIOD IN JULIAN DAYS 

/D1*D5/ 

B BERTH 

/I, 2/ 

10 INBOARD OR OUTBOARD BERTH 

/INBD, OUTBD, OUTBDl/ 

BASIN AREA BETWEEN TWO ADJACENT PIERS 

/BSl , BS2 , BS3 , BS4 , BS5 , BS6/ ; 

ALIAS (I,IP); 
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★ 

★ 

* 



GIVEN DATA 



SCALAR 

ADRAFT 

TUGWIDTH 

EXTEND 

BETWEEN 

FENDER 



SAFETY DISTANCE BET SHIP AND WATER BOTTOM /5/ 

BEAM WIDTH OF A STANDARD TUG /llO/ 

MAXIMUM DISTANCE SHIPS MAY EXTEND PAST PIER /20/ 

BOW STERN DISTANCE BETWEEN SHIPS AND OR SEAWALL /20/ 
FENDER DISTANCE BETWEEN NESTED SHIPS AND PIER /lO/; 



TABLE 



PIERDATA(P, PDATA) LIST OF PIERS AND CHARACTERISTICS 





PIERL 


DEPTH 


POWER 


12N 


1300 


50 


24 


12S 


1300 


50 


24 


IIN 


1397 


50 


24 


IIS 


1397 


50 


24 


ION 


1300 


38 


56 


7N 


1350 


45 


24 


7S 


1350 


45 


21 


4S 


1347 


40 


24 r 



TABLE 



SHIPDATA(I ,SDATA) 



LIST OF SHIPS AND CHARACTERISTICS 





LENGTH 


DRAFT 


BEAM 


ARRIVE 


DEPART 


SHIPPWR 


INSPECT 


DEPLOY 


SHP 


COSTSHF 


AFS2 


581 


24 


79 


5 


5 


4 








400 


AOR4 


659 


33 . 3 


96 


3 


5 


3 








500 


DDG4 


437 


20 


47 


2 


5 


3 








300 


DDG6 


437 


20 


47 


1 


5 


3 








300 


LHA4 


840 


26 


106 


2 


5 


14 






1 


400 


BB61 


887 


38 


108 


1 


2 


6 








500 


BB61a 


887 


38 


108 


5 


5 


6 








500 


CVN71 


1300 


37 


134 


1 


5 


8 






1 


1000 


CV67 


1300 


35.9 


130 


1 


5 


24 




200 


1 


1000 


CV66 


1300 


37 


130 


1 


3 


24 




200 


1 


1000 


CG27 


547 


28.8 


54 . 8 


1 


4 


4 








350 


CG30 


547 


28 . 8 


54 . 8 


2 


5 


4 








350 


CG34 


547 


28.8 


54 . 8 


4 


5 


4 








350 


CG4 8 


567 


33 


55 


1 


3 


6 


100 






350 


CG51 


566 


31 


55 


1 


2 


6 








350 


FFG5 


414 


24.2 


44 . 2 


1 


4 


2 




100 




300 


TAF8 


524 


22 


72 


1 


2 


2 








400 


TA0189 


678 


34.5 


97.5 


2 


4 


4 








4 00; 
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TABLE 



AVAIL(P, 


J) 






LIST OF PIERS 


WITH 


SERVICES OFFERED 




DFM 


JP5 


MSC 


STR 


140T 4160V DRON2 


ORDN 


CD68 TEND 


12N 


1 


1 






1 


1 




12S 


1 


1 






1 


1 




IIN 










1 


1 




IIS 










1 


1 




ION 














1 


7N 


1 


1 








1 


1 


7S 


1 










1 




4S 


1 


1 


1 


1 


1; 







TABLE 

LOGREQ(I,J) WEIGHTED SHIP TO SERVICES REQUIREMENTS 

DFM JP5 MSC STR 14 OT 4160V DR0N2 ORDN CD68 TEND 
AFS2 600 600 

A0R4 750 750 
LHA4 
DDG4 
DDG6 
BB61 

BB61a 600 



CVN71 


900 


999 




CV67 


900 






CV66 


900 






CG27 






400 


CG3 0 






400 


CG34 






400 


CG48 






400 


CG51 






400 



FFG5 200 

TAF8 700 700 400 

TA0189 700; 

TABLE 



BELONG (BASIN, P) PIERS BELONGING TO THE SAME BASIN 



BSl 

BS2 

BS3 

BS4 

BS5 

BS6 



12N 12S IIN IIS ION 7N 7S 4S 
1 

1 1 

1 1 

1 

1 

1 ; 
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TABLE 



DESIRE(I,P,B,IO,T) DESIRE POSITIONS FOR SHIPS BY DAY 



BB61. (7N,7S) . (1,2) . INBD 

CVN71. 12N. 2 . INBD 

CV66. 12S .1. inbd 

CV67 . (IIN, IIS) . (1,2) .INBD 

CG30.10N. (1,2) . INBD 

CG34 . ION. (1,2). inbd 

CG51.10N. (1,2) . (INBD,OUTBD) 



D1 D2 D3 D4 D5 
1 1111 
1 1111 
1 1111 
1 1111 
1 1111 
1 1111 
1 11 ; 



PARAMETER 

PREV(I,P,B,IO) PREVIOUSLY ASSIGNED SHIPS; 



PREV ( ' FFG5 ' , ' 1 IN ' , ' 2 ' , ' OUTBD ' ) =1 ; 
PREV( 'CVN71' , '12N' , '2 ' , 'INBD' ) =1 ; 
PREV( 'CG27 ' , ' IIN ' , ' 2 ' , ' INBD' ) =1; 
PREV( ' CG48 ' , ' ION ' , ' 2 ' , ' INBD' ) =1 ; 
PREV( 'CG51 • , ' ION INBD' ) =1 ; 
PREV( 'TAF8' , '4S' , '2' , 'INBD')=1; 
PREV ( ' CV66 ' , ' 12S ' , ' 1 ' , ' inbd ' ) =1 ; 
PREV ( ' DDG6 ' , ' 7S ' , ' 2 ' , ' INBD ' ) =1 ; 
PREV( 'BB61' , '7S' , '1' , 'INBD' ) =1 ; 



PARAMETER 

NOOUT(I) 

/CG27 1/ 

REWARD (10) 



IDENTIFIES SHIPS THAT CANNOT NEST SHIPS OUTBOARD 
REWARD FOR ASSIGNING SHIP INBOARD VICE OUTBOARD 



/INBD 

OUTBD 

OUTBDl 



300 

200 

100 / 



BWIDTH (BASIN) WIDTH OF BASIN 



/BSl 


365 


BS2 


765 


BS3 


700 


BS4 


283 


BS5 


365 


BS6 


2 2 0/ 
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FENSUP(I,P) FENDER AND SUPERSTRUCTURE RESTRICTION OR SUBJECTIVE LIMITS 



/ CVN71. ION 


-1 


CVN71.7N 


-1 


CV66, ION 




CV67, ION 


-1 


LHA4 ,45 


-1/; 



* 

* DERIVED DATA 

★ — 



SCALAR 

DAY COUNTER OF DAY /!/; 

PARAMETER 



INPORT(I) 

BENEFIT ( I, P,B, IO,T) 

CHKSPEC(I,T) 

COMPAT(I,P,B,IO,T) COMPATABLE SHIP-TO-PIER ASSIGNMENTS; 
SHIPDATA(I , 'SHP' ) = SHIPDATA ( I , ' SHP ' ) *2+1 ; 



CHKSPEC(I,T) = SUM( (P,B,IO), DESIRE ( I , P, B , 10 , T) ) ; 

COMPAT ( I , P , B , 10 , T) =0 ; 



COMPAT(I , P, B, IO,T) = 1 
$ ( ( SHIPDATA (I, 'DRAFT' ) 

AND ( SHI PDATA( I, 'LENGTH' ) 
AND (SHIPDATA (I, 'DEPART' ) 
AND ( SHIPDATA (I, 'ARRIVE' ) 
AND (FENSUP(I,P) 



LE (PIERDATA(P, 'DEPTH ')-ADRAFT) ) 
LE PIERDATA(P, ' PIERL' ) ) 

GE ORD(T)) 

LE ORD(T)) 

NE -1) ) ; 



COMPAT(I,P,B,IO,T)$(CHKSPEC(I,T) GT 0) = 1 
$( (DESIRE(I,P,B,IO,T) EQ 1) 

AND (COMPAT(I,P,B,IO,T) EQ 1)); 



COMPAT(I, P, B, 'OUTBD' ,T) $ (SHIPDATA(I , 'SHP' ) EQ 3)=0; 



COMPAT(I, P,B, 'OUTBDl ' ,T) $ (SHIPDATA(I , 'SHP' ) EQ 3)=0; 

BENEFIT (I,P,B,IO,T) $(COMPAT(I,P,B,IO,T) GE 1) = EXP ( - (ORD (T) -1 ) /7 ) * 

( SUM(J, LOGREQ(I,J) *AVAIL(P,J) ) 

+ SHIPDATA(I INSPECT' ) 

+ SHI PDATA( I, ' DEPLOY ' ) 

+ REWARD(IO)); 



INPORT(I)=SUM( (P,B,IO) , PREV (I , P, B, 10) ) ; 
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* 

DECISION VARIABLES 

* 



BINARY VARIABLES S ( I , P , B, lO , T) ; 

POSITIVE VARIABLES 
ZP(I,P,B,IO,T) ; 

VARIABLES UTILITY; 



FORJIULATION 



EQUATIONS 

SUCCESS EQUATION THAT MEASURES THE SUCCESS OF BERTHING SCHEDULE 
PIERLEN(P,T) SUM OF SHIPS LENGTH MUST BE LESS THAN PIER LENGTH 
POWER(P,T) PIER POWER CABLE EQUATION FOR EACH TIME PERIOD 
SHIPCOM(P, B,T) OUTBOARD SHIP LENGTH MUST BE LESS INBOARD SHIP 
SHIPCOMl (P, B,T) OUTBOARDl SHIP LENGTH LESS RHAN OUTBOARD SHIP 
ENSURE(I,T) ASSIGN EACH SHIP TO ONE BERTH IF THE SHIP IS IN 
ENSUREl (P, B,T) ASSIGN AT MOST THREE SHIPS PER BERTH 
ENSURE2 (P, B, IO,T) ASSIGN ONLY ONE SHIP TO ONE SPOT 
CHKDAYl(I,P,B,IO) PENALTY FOR BERTH SHIFT FOR DAY 1 
CHECKl(I,P,B,IO,T) PENALTY FOR BERTH SHIFT 

OUTBLIM(I,P,B,IO,T) ENSURES IDENTIFIED SHIPS HAVE NOONE OUTBOARD 
BASINLIM(BASIN, B,T) WIDTH OF SHIPS PLUS ROOM FOR TUG IN BASIN; 

SUCCESS. . 

UTILITY =E= SUM((I,P,B,IO,T)$COMPAT(I,P,B,IO,T) , 

S(I,P,B,IO,T) * BENEFIT ( I, P, B, IO,T) ) 

- SUM( (I , P, B, IO,T) $ ( (COMPAT(I, P, B, IO,T) EQ 1) AND 
(SHIPDATA(I, 'ARRIVE' ) LT ORD(T))), 

SHIPDATA(I, 'COSTSHF' ) *ZP(I , P, B, IO,T) ) 
-SUM((I,P,B,IO)$(COMPAT(I,P,B,IO, 'Dl') EQ 1 AND 
PREV(I , P, B, 10) EQ 1) ,SHIPDATA(I, 'COSTSHF' ) *ZP ( I , P, B, lO , ' Dl ' ) ) ; 



PIERLEN(P,T) . . 

SUM( (I , B) $COMPAT( I, P, B, ' INBD' ,T) , 

( SHIPDATA(I, 'LENGTH' )+BETWEEN) *S(I,P,B, 'INBD' ,T) )-Sl(P,T) 
=L= PIERDATA(P, 'PIERL')+EXTEND; 



CHKDAY1(I,P,B,I0)$(C0MPAT(I,P,B,I0, 'Dl') GE 1 
S(I,P,B,IO, 'D1')+ZP(I,P,B,I0, 'Dl') =G= 



AND PREV(I,P,B,IO) 

EQ 



PREV(I,P,B,IO) ; 



1 ) . . 
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CHECKl(I,P,B,IO,T)$(ORD(T) GT 1 AHD 

SHIPDATA (I, 'ARRIVE' ) LTORD(T) AND 

COMPAT(I, P, B, IO,T) GT 0).. 



S(I,P,B,IO,T)-(COMPAT(I,P,B,IO,T-l)*S(I,P,B,IO,T-l) ) -ZP ( I , P, B , lO, T) 

=L=0; 



POWER (P,T) . . 



SUM( (I,B,IO)$COMPAT(I,P,B,IO,T) , 

S (I,P, B, IO,T) *SHIPDATA(I, 'SHIPPWR* ) ) -S2 (P, T) =L=PIERDATA (P, 'POWER' ) 
SHIPCOM(P,B,T) 

$( SUM(I,COMPAT(I,P,B, 'OUTBD' ,T) ) GT 0).. 

SUM(I$COMPAT(I,P,B, 'INBD' ,T) , 

S (I,P,B, 'INBD' ,T) *SHIPDATA(I, 'LENGTH' ) ) =G= 

SUM(I$COMPAT(I,P,B, 'OUTBD', T) , 

S(I,P,B, 'OUTBD' ,T) *SHIPDATA(I, 'LENGTH') ) ; 



SHIPC0M1(P,B,T) 

$(SUM(I,COMPAT(I,P,B, 'OUTBDl' ,T) ) GT 0).. 

SUM(I$COMPAT(I,P,B, 'OUTBD' ,T) , 

S (I, P,B, 'OUTBD' ,T) *SHIPDATA(I, 'LENGTH ' ) ) =G= 
SUM(I$COMPAT(I,P,B, 'OUTBDl' ,T) , 

S (I,P,B, 'OUTBDl ' ,T) *SHIPDATA(I, 'LENGTH' ) ) ; 



ENSURE(I,T)$(SUM( (P,B,IO) ,COMPAT(I,P,B,IO,T) ) GE 1).. 

SUM((P,B,IO)$COMPAT(I,P,B,IO,T) ,S(I,P,B,IO,T)) =E= 1; 



ENSUREl (P, B,T) 

$(SUM((I,IO) ,COMPAT(I,P,B,IO,T)*(SHIPDATA(I, 'SHP')-!)) GT 2).. 

SUM( (I,IO)$COMPAT(I,P,B,IO,T) , 

SHIPDATA(I, 'SHP') *S(I,P,B,IO,T) ) =L= 3; 

ENSURE2(P,B,IO,T)$(SUM(I,COMPAT(I,P,B,IO,T)) GT 1).. 
SUM(I$COMPAT(I,P,B,IO,T) , S (I , P, B, IO,T) ) -S5 (P, B, IO,T) =L= 1; 
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BASINLIM ( BASIN, B,T)$ (SUM (P, BELONG ( BASIM, P) ) GT 0 
AND SUM( (I,P,IO) ,COMPAT(I,P,B,IO,T)*BELONG(BASIN,P) ) GT 1)., 

SUM( (I,P,IO) $COMPAT(I,P,B,IO,T) , 

S(I , P, B, IO,T) * (SHIPDATA(I, 'BEAM' ) +FENDER) * BELONG ( BASIN , P) ) 

=L= BWIDTH(BASIN)-TUGWIDTH; 

OUTBLIM(I,P,B,IO,T) $(NOOUT(I) EQ 1 AND ORD(IO) LT 3 
AND COMPAT(I , P, B, IO,T) EQ 1).. 

(3-ORD(IO) ) *S (I , P, B, IO,T) -COMPAT ( I , P , B , lO+l , T) *S (I , P, B, IO+l,T) 
-COMPAT ( I , P, B , 10+2 , T) *S ( I , P, B , 10+2 , T) 

+SUM(IP $COMPAT(IP, P, B, I0+1,T) , S ( IP , P, B, lO+l , T) ) 

+SUM(IP $C0MPAT(IP,P,B,I0+2,T) , S ( 1 P, P , B , 10+2 , T) ) =L= ( 3-ORD ( lO) ) 



MODEL SHIPl /ALL/; 

OPTION LIMROW=0, LIMCOL=0, OPTCR=.0, ITERLIM=50000 , RESLIM=2500; 
OPTION SOLPRINT = On, SYSOUT = Off; 

SOLVE SHIPl USING MIP MAXIMIZING UTILITY; 

DISPLAY S.L; 

PARAMETER SOL(I , P, B, 10) ; 

OPTION SOL: 1:1: 3; 

DAY = 1; 

LOOP(T, 

SOL(I,P,B,IO) = S.L(I,P,B,IO,T) ; 

DISPLAY DAY, SOL; 

DAY = DAY +1) ; 
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